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摘要 :气候 变化 已 明显 影响 到 陆地 植被 的 活动 ,但 在 不 同 生态 系统 间 存 在 差异 ,人 研究 不 同 陆 地 生态 系统 归 一 化 植被 指数 ( NDVI) 
的 时 空 变化 特征 ,不 仅 可 揭示 各 生态 系统 植被 活动 对 气候 变化 的 响应 规律 ,而 且 可 为 我 国 不 同 生 态 区 制定 应 对 气候 变化 的 策略 
和 生态 文明 建设 提供 科学 依据 。 基 于 1982 一 2012 年 CIMMS NDVI3g 和 中 国 陆 地 生态 系统 类 型 数据 ;利用 一 元 线性 回归 集合 
经 验 模 态 分 解 和 相关 分 析 等 方法 ,研究 了 近 30 年 中 国 各 陆地 生态 系统 NDVI 的 时 空 变化 特征 ,分 析 了 其 与 气候 事件 的 关系 。 
结果 表明 , 近 30 年 中 国 植被 活动 显著 上 升 ,年 平均 归 一 化 植被 指数 (ANDVI) 的 睫 升 幅度 为 0.0029/10a( P<0.05) ,年 最 大 归 一 化 
植被 指数 (MNDVI) 的 上 升幅 度 为 0.0076/10a( P<0.01) ;植被 活动 显著 增强 的 区 域 主要 是 分 布 在 东部 季风 区 的 农田 和 森林 生态 
系统 ,显著 下 降 的 区 域 主要 是 分 布 于 西北 的 蕊 漠 生 态 系 统 和 东北 的 森林 生态 系统 ;尽管 ANDVI 和 MNDVI 线性 趋势 的 显著 性 有 
所 差异 ,但 农田 森林 、 草 地 和 水 体 与 湿地 生态 系统 的 NDVI 总 体 呈 非 稳定 的 上 升 趋势 ,上 升 过 程 中 伴随 着 较 大 波动 , 藉 漠 生态 
系统 的 NDVI 呈 下 降 趋 势 , 植 被 退化 显著 ; 与 线性 趋势 不 同 ,各 生态 系统 植被 活动 的 残 差 趋势 包含 “上升 一 下 降 ”两 个 阶段 ,并 相 
继 于 20 世纪 90 年 代 到 21 世纪 初 发 生 转折 ;上 述 5 类 生态 系统 的 植被 活动 存在 不 同 尺 度 的 周期 特征 ,年 际 周期 波动 特征 (1.9 一 
7.6a) 比较 显著 ,而 年 代 际 周期 (10.7a 和 22.2a) 的 显著 性 相对 较 差 ;各 生态 系统 的 空间 异 质 性 在 趋 强 过 程 中 存在 2.1 一 7.1a 的 年 
际 周期 节律 ; 海 详 与 大 气 环流 的 短 周期 脉动 与 各 生态 系统 植被 活动 的 周期 性 节律 有 着 明显 关联 ,ENSO 事件 和 太阳 活动 是 推动 
植被 活动 周期 性 振荡 的 重要 因素 。 
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Abstract: Vegetation, an important component of terrestrial ecosystems, tells the story of the ecological response to climate 
change. In-the context of climate change, vegetation changes have directly affected substance exchanges and energy balance 
of regional land-atmosphere interaction processes. Additionally, in China, the vegetation activity of varied terrestrial 
ecosystems has responded to global climate change in significantly different ways. As a vital remote sensing indicator, the 
normalized difference vegetation index ( NDVI) reflects vegetation growth and coverage, and time series analysis can 


describe the temporal evolution and spatial variability of vegetation very well. Thus, studying the spatiotemporal variation of 
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NDVI in terrestrial ecosystems not only reveals the response of vegetation in those ecosystems to climate change, but can 
provide a scientific basis for action on climate change and ecological considerations in an increasingly urbanized country like 
China. Based on the Global Inventory Modeling and Mapping Studies (GIMMS) NDVI3g dataset for 1982—2012, and the 
macro structure data of the Chinese terrestrial ecosystem, this study utilized one-dimensional linear regression, ensemble 
empirical mode decomposition (上 EEMD ) ，and correlation analysis to examine spatiotemporal variation in the NDVI of each 
terrestrial ecosystem in China. lt also analyzed the relationship between NDVI and climate events over a recent 30:year 
period. The results show that the average normalized difference _ vegetation index ( ANDVI ) and maximum normalized 
difference vegetation index ( MNDVI) increased (0.0029/10 a, P<0.05 and 0.0076/10 a, P<0.01, respectively ) over the 
same period. These phenomena indicate the activity of vegetation strengthened in China. In addition, owing to its sensitivity 
to water-thermal factors, MNDVI can be a useful supplement for monitoring vegetation change. Furtheimiore, the areas with 
significant increases in vegetation activities were mainly the farmland and forest ecosystems in the eastern molisoon region, 
and the areas with significant decreases in vegetation activities were mainly the forests in Northeast China ‘and the desert 
ecosystem in Northwest China. With ANDVI and MNDVI showing linear trends at different significance levels, the NDVIls of 
the farmland, forest, grassland, and water and wetland ecosystems showed a general upward trerid with large fluctuations, 
while the NDVI of the desert ecosystem decreased, with obvious vegetation degradation However, EEMD detected an 
inflection point among the trends in the 和 ve ecosystems that was distinct from'thé"linear trend. This occurred in the 1990s 
and early 21” century when the NDVI trend shifted from increasing to decreasing The decomposed results indicated that the 
NDYVI of the ecosystems had oscillation characteristics on an inter-annual (1:9—7.6years) and an inter-decadal (10.7 一 
22.2 years ) scales, with the former being statistically more prominent、The increasing spatial heterogeneity of the NDVI of 
different terrestrial ecosystems also showed inter-annual periods of 2.1 一 2. /years ，with increasing regional differences in 
vegetation activity. Furthermore, the short period-pulsatiomof the'alirssea system from the El Nifo-Southern Oscillation and 
solar activity influenced vegetation activity via climatic factors ，whieh was evidently related to the oscillation characteristics 


in all ecosystems. 
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气候 变化 已 影响 到 陆地 生态 系统 的 生物 物理 和 化 学 过 程 "” ,植被 作为 陆地 生态 系统 的 重要 组 成 部 分 ， 
对 气候 变化 比较 敏感 ,气候 变化 引起 的 植被 变化 已 直接 影响 到 区 域 陆 一 气 交 互 过 程 的 物质 与 能 量 平 
衡 '" 。 因 此 ,开展 区 域 及 全 球 尺度 的 植被 变化 研究 对 认识 气候 变化 与 陆地 生态 系统 的 相互 作用 有 重要 意 
义 “。 大 尺度 植被 变化 的 连续 观测 对 传统 生态 学 研究 技术 提出 了 挑战 , 随 着 全 球 变 化 生态 学 的 发 展 ,遥感 
技术 逐渐 成 为 大 尺度 生态 问题 连续 观测 的 重要 手段 。 卫 星 传感器 通过 探测 陆 表 的 光谱 反射 信息 来 识别 植被 
的 变化 ,利用 植被 在 可 见 光 与 近 红 外 波段 的 反射 特征 构建 的 归 一 化 植被 指数 (Normalized Difference Vegetation 
Index, NDVI) , 常 作为 植被 生长 状况 及 覆盖 度 变 化 的 最 佳 指 示 因 子 ” ,已 在 不 同 尺 度 植被 变化 研究 中 广泛 应 
用 “*" ,特别 是 时 间 序列 的 NDVI 能 较 好 的 反映 陆 表 植被 活动 的 时 间 演 化 和 空间 变异 特征 。Myneni 等 
利用 1982 一 1991 年 的 AVHRR-NDVI 数据 对 全 球 植被 活动 进行 了 研究 ,提出 了 北半球 中 高 纬 地 区 植被 活动 增 
强 的 观点 ;Piao 等 "中 利用 Global Inventory Modeling and Mapping Studies (GIMMS) 第 二 代 NDVI 产品 研究 了 
1982 一 2006 年 亚 欧 大 陆 生长 季 的 植被 活动 ,指出 NDVI 存在 两 个 明显 的 变化 阶段 ,1997 之 前 显著 上 升 ,1997 
年 之 后 NDVI 上 升 趋势 减缓 ;Park 等 ' 外 同样 利用 该 CIMMS 产品 研究 了 东亚 地 区 的 植被 活动 ,得 出 NDVI 增加 
趋势 在 20 世纪 90 年 代 中 期 发 生 明 显 转 折 ; 杜 家 强 等 利用 AVHRR 和 MODIS 两 个 传感器 构建 的 连续 NDVI 
数据 对 新 疆 近 30 年 的 植被 活动 展开 研究 ,认为 NDVI 的 变化 趋势 在 20 世纪 90 年 代 末 存在 拐点 。 此 外 ,诸多 
学 者 关于 不 同 地 理 单元 …… 力 至 全 国 范围 ”的 植被 活动 研究 认为 ,人 研究 区 的 植被 活动 总 体 增强 。 
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然而 ,以 上 研究 存在 几 方 面 的 问题 ,一 是 大 多 研究 基于 上 一 代 GIMMS NDVIL2( 1982 一 2006) 、SPOT-VGT 
NDVI(1998 一 ) 或 MODIS NDVI(2000 一 至 今 ) 等 数据 ,往往 研究 时 间 序 列 较 短 ;尽管 有 学 者 延续 不 同 传感器 的 
NDVI 数据 开展 植被 变化 研究 ,但 不 同 传感器 的 性 能 差异 又 增加 了 NDVI 数据 的 不 确定 性 。 二 是 传统 的 线 
性 分 析 方 法 对 时 序 较 短 的 非 平稳 NDVI 序 列 开展 变化 趋势 研究 ,会 忽略 非 平稳 数据 的 非 线性 波动 特征 ,致使 
人 研究 结论 和 规律 存在 局 限 性 ,而 集合 经 验 模 态 分 解 ( Ensemble Empirical Mode Decomposition , EEMD ) 因 其 良好 
的 自 适应 性 可 弥补 该 缺陷 。 三 是 前 人 研究 多 从 自然 单元 和 行政 区 划 的 角度 探讨 植被 活动 特征 ,忽视 了 不 同 陆 
地 生态 系统 内 部 的 植被 活动 规律 ,而 我 国 类 型 多 样 的 陆地 生态 系统 往往 灸 拣 分 布 , 不 同 陆地 生态 系统 的 植被 
活动 对 气候 变化 的 响应 存在 明显 差异 。 为 此 ,本 研究 将 利用 1981 一 2012 年 的 第 三 代 GIMMS NDVI 数据 
(GIMMS NDVI3g) 和 集合 经 验 模 态 分 解 等 非 线 性 分 析 方法 ,从 生态 系统 的 视角 出 发 ,开展 近 30 年 中 国 陆 地 生 
态 系统 NDVI 时 空 变化 特征 研究 ,以 期 揭示 气候 变化 背景 下 各 生态 系统 的 啊 应 ,为 我 国 制定 应 对 气候 变化 的 
策略 和 生态 文明 建设 提供 科学 依据 。 


1 数据 与 方法 


1.1 数据 资料 
1.1.1 GIMMS NDVI3g 数据 

美国 国家 航空 航天 局 (National Aeronautics and Space Administration, NASA) 基 于 NOAA-AVHRR 系列 数据 
开发 出 的 GIMMS NDVI3g 产品 是 一 种 长 时 序 的 全 球 植被 指数 数据 集 “(1981 二 2012 年 ) ,空间 分 辨 率 为 8km， 
时 间 分 辨 紊 是 13d。 该 数据 集 消除 了 大 气 水 汽 、 火山 爆发 ,太阳 高 度 角 和 传 感 希 灵敏 度 的 变化 等 因素 的 影 
响 “ ,精度 得 到 有 效 提高 ,良好 的 数据 质量 及 连续 的 时 间 序 列 使 其 在 大 尺度 生态 过 程 模拟 研究 中 得 到 广泛 应 
用 。 本 人 研究 从 NASA 网 站 (http://ecocast.arc.nasa.gov/ ) 获取 了 1982 一 2012 年 覆盖 全 球 的 GIMMS NDVI3g 数 
据 , 运 用 IDL 软件 编程 对 膛 感 影像 进行 读 取 、 配 准 、 格 式 转 换 以 及 研究 区 裁剪 等 预 处 理 。 
1.1.2 中 国 陆地 生态 系统 类 型 数据 

中 国 陆地 生态 系统 类 型 宏观 结构 数据 源 于 中 国 科学 院 资 源 环 境 科 学 数据 中 心 ( http://www.resdc.cn/)。 
基于 遥感 解 译 获取 的 1:10 万 土地 利用 /十 地 覆盖 数据 对 各 生态 系统 类 型 进行 辨识 研究 和 重新 分 类 ,形成 中 
国 陆 地 生态 系统 类 型 数据 集 ( 分 别 为 农田 生态 系统 、 和 森林 生态 系统 .草地 生态 系统 ,水体 与 湿地 生态 系统 .项 
澳 生 态 系统 .聚落 生态 系统 和 其 他 生态 系统 ) 。 本 研究 采用 2010 年 的 数据 集 ,提取 不 同 生态 系 统 类 型 , 因 聚 落 
生态 系统 和 其 他 生态 系统 类 型 所 占 比 例 较 低 ,同时 考虑 到 聚落 生态 系统 的 碎 小 斑 块 很 难 在 8km 分 辩 率 的 
GIMMS NDVI3g 数据 中 反映 出 来 ,其 他 生态 系统 主要 为 裸 土地 或 裸 宕 石 质地 ,其 NDVI 指标 无 实际 生态 学 意 
义 ,故人 研究 中 未 对 这 两 类 生态 系统 类 型 进行 分 析 。 
1.1.3 其 他 数据 

多 变量 厄尔尼诺 -南方 涛 动 指数 (The Multivariate ENSO Index，MEI) 来 源 于 美国 国家 海洋 与 大 气管 理 局 
地 球 系 统 人 研究 实验 室 (http://www.cdc.noaa.gov/people/klaus.wolter/ MEI/) ,集合 了 海平 面 气压 海 表 温度 、 
面 气温 地面 第 向 风 和 径 向 风 、 总 云 量 等 参数 ,能 够 全 面 地 监测 和 反映 厄尔尼诺 -南方 涛 动 ( El Nifio-Southern 
Oscillation ，ENSO ) 事件 ””。 太 阳 黑 子 数据 采用 比利时 皇家 天 文 台 太阳 活动 数据 分 析 中 心 ( http://www.sidc. 
be/silso/datafiles ) 发 布 的 月 太阳 黑子 数 序列 ,该 数据 已 采用 13 个 月 滑动 窗口 平滑 处 理 。 本 研究 截取 1982 一 
2012 年 的 时 间 序 列 , 并 用 集合 经 验 模 态 分 解 提取 不 同时 间 尺 度 的 周期 分 量 。 
1.2 研究 方法 
1.2.1 年 植被 指数 合成 方法 

中 国 植被 类 型 丰富 多 样 ,同一 生态 系统 往往 兼 有 多 个 目 然 齐 的 不 同 植被 型 , 且 部 分 植被 型 存在 明显 的 季 
相 特征 。 为 更 全 面 描述 年 际 植被 活动 特征 ,本 研究 合成 了 年 平均 归 一 化 植被 指数 (ANDVI) 和 年 最 大 归 一 化 
植被 指数 ( MNDVI) 两 种 指标 ” ,分 析 1982 一 2012 年 中 国 陆地 生态 系统 植被 活动 的 空间 变化 。 其 中 ,ANDVI 
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反映 年 内 植被 生长 的 平均 状态 ,MNDVI 反映 当年 最 优 水 热 组 合 下 的 植被 最 佳 生长 状态 , 即 最 佳 植 被 长 势 。 


1 24 

ANDVI = 27 > NDVI(1 < i < 24,1982 <j < 2012) (1) 
i=24 

MNDVI = max( NDVL) (1 < i < 24,1982 </ < 2012) (2) 


1.2.2 ”回归 分 析 

采用 最 小 二 乘法 逐 像 元 拟 合 1981 一 2012 年 间 的 CIMMS NDVI3g 数据 序列 , 求 取 各 像 元 近 30 年 的 变化 斜 
率 。 一 元 线性 回归 和 斜率 值 的 正 、 负 分 别 表示 NDVI 升 高 或 降低 的 趋势 ,斜率 的 大 小 则 反映 了 NDVI 值 工 升 或 下 
降 的 程度 ,回归 趋势 的 显著 性 采用 下 检验 判定 。 将 NDVI 变化 趋势 和 不 同 显著 性 水 平 (P=0.05 和 P=0.01) 羞 
加 分 析 , 可 得 到 极 显 著 上 升 . 显 著 上 升 上 升 (但 不 显著 ) .下 降 ( 但 不 显著 ) .显著 性 下 降 和 极 显 著 性 下 降 6 种 
变化 趋势 组 合 。 
1.2.3 集合 经 验 模 态 分 解 

集合 经 验 模 态 分 解 (Ensemble Empirical Mode Decomposition, EEMD) 2 是 近年 来 发 展 的 一 种 适应 于 非 
线性 、 非 平稳 序列 的 信号 分 析 方 法 ,该 法 在 EMD( Empirical Mode Decompoesition ) 的 基础 上 ,引入 辅助 白 噪声 ， 
继承 了 EMD 的 上 自 适 性 ,有 效 地 解决 了 尺度 混合 问题 ; 它 将 原始 信号 中 不 同时 间 太 度 的 周期 特征 分 解 开 , 形 成 
具有 不 同 尺 度 周 期 的 本 征 模 态 函 数 分 量 ( Intrinsic Mode Function，IMF) 和 一 个 残 差 分 量 ( Residual Standard 
Error，RSE) 。 其 中 ,所 添加 白 噪 声 对 原始 信号 的 影响 将 遵循 如 下 统计 规律 : 


rc=e/ VN 或 Ing + 六 nN (3) 
式 中 , o 为 标准 偏差 , 即 输 入 信号 与 相应 IMF 分 量 的 偏差 ,s 为 自 噪声 幅度 , V 为 集合 样本 数 。 
经 蒙特 卡 罗 法 验证 , 当 自 由 度 趋 于 无 穷 ,添加 白 噪 声 的 序列 经 过 分 解 后 ,IMF 的 能 量 谱 平 均 密度 与 平均 周 
期 具有 如 下 关系 : 


lInk, + ln7, =0 (4) 
式 中 , ,为 n 个 IMF 的 能 量 谱 平均 密度 ,7 为 相应 平均 周期 。 理 论 上 , 白 噪声 IMF 分 量 的 平均 周期 与 谱 密 度 
在 二 维 空间 中 分 布 在 斜率 为 -1 的 直线 上 ,但 实际 偏差 使 得 连续 概率 分 布下 日 噪声 能 量 详 密 度 与 平均 周期 有 


如 下 关系 : 
Fa 7 2 ln7,/2 
Ink, = ln7T, +k a 5 (5) 


式 中 , ,为 平均 能 量 ,了 了, 为 平均 周期 ,df 为 自由 度 , 当 有 = 2.326 和 1.645 分 别 可 以 得 到 99% 和 95% 的 置信 区 
间 。 分 布 在 不 同 显著 水 平 曲线 之 上 的 IMF 分 量 表示 通过 了 相应 水 平 的 显著 性 检验 , 即 该 IMF 分 量具 有 显著 
的 物理 意义 。 

EEMD' 分 解 15d 合成 的 GIMMS NDVI3g 序列 及 标准 差 序列 时 ,辅助 白 品 声 与 原始 信号 的 标准 差 比 率 为 
0.2 ,集合 样本 数 为 1000 ,并 借助 白 噪 声 的 集合 扰动 对 IMF 分 量 进行 显著 性 检验 ,获取 不 同 IMF 的 信和 度 水 平 。 
其 中 ,IMF 分 量 表 征 从 高 频 到 低频 不 同时 间 尺 度 的 NDVI 周期 波动 特征 ,RSE 反映 了 NDVI 时 间 序 列 的 残余 变 
化 趋势 。 


2 ， 结果 与 分 析 


2.1 近 30 年 中 国 NDVI 总 体 空间 变化 特征 

1982 一 2012 年 ,中国 ANDVI 整体 旺 上 升 趋势 (0.0029/10a,P<0.05)。 从 图 1 看 出 ,全 国有 59.31% 的 国土 
区 域 ANDVI 呈 上 升 趋势 ,其 中 仅 44.69% 的 区 域 超过 5% 显著 性 水 平 ;全 国有 40.69% 的 国土 区 域 ANDVI 呈 下 
降 趋势 ,其 中 有 46.43% 的 区 域 超过 5% 显著 性 水 平 。 从 空间 上 看 ,大 兴安 岭 .嫩江 流域 上 游 、 三 江平 原 、 长 白 
山 .长 三 角 地 区 和 珠 三 角 地 区 的 ANDVI 显著 下 降 , 降 幅 达 到 0.0080/10a 以 上 (P<0.01) ,说 明 这 些 区 域 的 年 平 
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均 植 被 活动 在 减弱 ;此 外 ,西北 地 区 的 准噶尔 盆地 .塔里木 盆地 . 柴 达 木 盆地 .甘肃 西部 和 内 蒙古 西部 广大 地 区 
的 ANDVI 降幅 不 大 ,但 下 降 趋 势 极 为 明显 ,这 与 上 述 地 区 NDVI 本 底 值 较 低 有 关 。 华 北平 原 .黄土 高 原 . 秦 巴 
山区 台湾 山脉 .东南 丘陵 局 部 ,海南 岛 局 部 和 横断 山区 的 ANDVI 极 显著 上 升 ,增幅 达到 0.0080/10a 以 上 (P< 
0.01) ,局 部 地 区 增幅 远 超 0.0240X10a ,说 明 这 些 地 区 的 年 平均 植被 活动 在 增强 ;此 外 , 藏 北 高 原 及 天 山 局 部 地 
区 的 ANDVI 虽然 增幅 不 高 ,但 上 升 趋势 极为 明显 。 


ANDVI 趋 势 /10a | / ANDVLF 检 验 » / 
国 <-0.008 | 国 极 显著 下 降 x a 
国 -0.008 一 0 | 国 显著 下 降 
国 0 一 0.008 1 “下 降 但 不 显著 a 
国 0.008 一 0.016 i 一 i 
0 Su 


0 1000 km 

MNDVI 趋 热 /10 a | MNDVI-F 检 验 上 / 
国 <-0.008 人 国 极 显著 下 降 、 0 
国 -0.008 一 0 国 显著 下 降 

国 ] 0 一 0.008 | “下 降 但 不 显著 EE i | 
国 | 0.008 一 0.016 , se 一 _ Ea 
国 0.016 一 0.024 攻 

男 | > 0.024 国 极 显 著 上 升 \ 4 


图 1 最 大 归 一 化 植被 指数 (ANDVI) .平均 归 一 化 植被 指数 ( MNDYVI) 变化 趋势 及 显著 性 检验 
Fig.l Trend and significance test of ANDVI and MDNVYI 


近 30 年 全 国 MNDVI 整 体 也 呈 上 升 趋势 ,上 且 上 升幅 度 比 ANDVI 更 强 (0.0076/10a,P<0.01)。 如 图 1 所 
示 , 占 比 65.24% 的 国土 区 域 MNDVI 呈 增 加 趋势 ,其 中 46.85$% 的 区 域 超过 $% 显著 性 水 平 ; 占 比 34.76% 的 国 
土 区域 MNDVI 呈 下 降 趋势 ,其 中 超过 5$% 显 著 性 水 平 的 区 域 仅 占 39.27%。 从 空间 上 看 , 极 显著 下 降 的 区 域 主 
要 分 布 在 以 荒漠 生态 系统 为 代表 的 内 蒙古 高 原 西部 .准噶尔 盆地 和 塔里木 盆地 ,说明 这 些 地 区 年 最 佳 植被 长 
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势 在 逐年 减弱 ;而 东北 平原 .华北 平原 .黄土 高 原 .天 山南 北 两 侧 绿洲 区 .大巴 山脉 两 广 丘陵 及 云 贯 高 原 局 部 
的 MNDVI 呈 极 显 著 上 升 趋势 ,增幅 达到 0.0240/10a 以 上 (P<0.01) ,说 明 这 些 地 区 年 最 佳 植被 长 势 在 逐年 
增强 。 

由 于 MNDVI 和 ANDYVI 代表 的 植被 活动 意义 不 同 , 导 致 两 者 的 变化 幅度 差异 较 大 ,尤其 在 东部 季风 区 ， 
MNDYVI 的 变化 幅度 明显 高 于 ANDVI。 同 时 ,二 者 的 变化 趋势 在 一 些 区 域 存在 差异 ,其 至 截然 相反 ,其 中 东南 
沿海 的 浙 六 丘陵、 台湾 名 和 海南 岛 地 区 的 年 平均 植被 活动 在 增强 ,但 年 最 佳 植 被 长 势 却 在 变 差 ,人 类 活动 持续 
增强 导致 森林 植被 破坏 可 能 是 重要 原因 ; 东北 松 嫩 平原 的 年 平均 植被 活动 在 减弱 ,这 与 东北 地 区 农田 
NDVI 下 降 趋 势 一 致 , 但 年 最 佳 植 被 长 势 却 在 变 好 ,这 种 分 异 可 能 是 暧 干 化 趋势 下 松 嫩 平原 农田 生态 系统 的 
种 植 结构 和 作物 熟 制 变 化 引起 的 。 结 合 ANDVI 和 MNDVI 来 看 ,全 国 植被 活动 显著 增强 的 区 域 主要 是 分 布 在 
东部 季风 区 的 农田 和 森林 生态 系统 ,而 西北 地 区 的 充 浇 生态 系统 和 东北 森林 区 的 植被 活动 则 显 若 减弱 
2.2 不 同 陆 地 生态 系统 的 NDVI 线性 变化 趋势 分 析 

中 国 各 生态 系统 面积 占 比 由 高 到 低 依 次 是 草地 森林、 农田、 充 漠 和 水 体 与 湿地 生态 系统 ; ANDVI 和 
MNDVI 最 高 的 是 森林 生态 系统 ,分 别 为 0.5297 和 0.8437 ,荒漠 生态 系统 最 低 ,分别 为 0.0789 和 0.1543。 线 性 
趋势 ( 表 1) 表 明 , 除 范 漠 生态 系统 外 ,各 生态 系统 的 植被 活动 主要 表现 出 上 升 趋势 农田 生态 系统 的 植被 活 
动 显 著 增强 ,ANDVI 的 增幅 为 0.0061/10a( P<0.01) ,MNDVI 达到 0.0162/10a( P<0.05) ,这 与 模型 模拟 的 中 
农业 植被 净 初 级 生产 力 增加 趋势 一 致 , 近 30 年 中 国 农业 种 植 结构 调整 \ 化 肥 用 量 增加 、 大 气 CO, 浓度 升 高 
和 气候 变 暧 “I 。 森 林 生 态 系统 的 ANDVI .MNDVI 均 呈 不 同 幅度 的 十 升 趋势 ;但 ANDVI 表现 出 了 较 强 的 空间 
差异 ,东北 地 区 和 森林 生态 系统 的 ANDVI 显著 下 降 ,而 南方 亚热带 森林 生态 系统 的 ANDVI 则 显著 上 升 ,这 与 森 
林 生 态 系 统 的 碳 源 / 汇 空间 分 布 格局 较为 一 致 ”1 ,MNDVI 极 显著 十 升 ,年 最 好 植被 长 势 显 著 增强 可 能 与 气候 
变化 有 关 ,气温 和 大 气 CO, 浓度 的 升 高 促进 了 和 森林 植被 年 内 最 成 盛 时 期 的 长 势 。 草 地 生态 系统 的 ANDVI 显 
著 上 升 (0.0022/10a,P<0.05) ,MNDVI 极 显 著 上 上 升 (0.0067/10a,P<0.01) ,说 明 我 国 草 地 地 上 生物 量 显著 增 
加 ,这 与 沈 海 花 等 “的 研究 结果 一 致 ,尽管 草地 生物 量 在 近 30 年 的 变化 趋势 存在 两 个 拐点 ,但 总 体 依然 呈 增 
加 趋势 。 水 体 与 湿地 生态 系统 的 NDVI 值 是 二 者 的 混合 信息 , NDVI 变化 主要 源 于 湿地 植被 变化 ,由 于 近 30 
年 我 国 湿 地 大 面积 萎缩 ,特别 是 水 体面 积 的 萎缩 和 湿地 植物 由 水 生 向 陆 生 的 演 替 , 均 会 造成 该 生态 系统 植 
被 盖 度 和 生物 量 的 增加 ,导致 NDVL 升 高 3 蒲 漠 生态 系统 与 前 4 类 生态 系统 截然 相反 ,ANDVI 和 MNDVI 都 处 
于 下 降 趋势 ,这 与 中 国 西北 日 益 严 重 的 落 漠 化 现状 比较 吻合 。 


表 1 各 陆地 生态 系统 近 30 年 ANDVI、MNDVI 基本 特征 及 变化 趋势 
Table 1 Basic statistics of ANDVI and MNDYVI in different terrestrial ecosystems and their trends in recent 30 years 


平均 归 一 化 植被 指数 ( ANDVI) 最 大 归 一 化 植被 指数 ( MNDVI) 
生态 系统 类 型 面积 比例 变化 幅度 变化 幅度 
Types of ecosystem Area 均值 方差 Amplitude 趋势 显著 性 均值 方差 Amplitude 趋势 显著 性 
2 proportion/ % Average Variance ofvariation/ Trends Significant Average Variance ofvariation/ Trends Significant 
(/10a) (/10a) 
农田 Farmland 18.81 0.4177 0.0163 0.0061 瞪 逢 六 0.7473 0.0188 0.0162 i 
和 森林 Forest 23.64 0.5297 0.0168 0.0053 上 升 0.8437 0.0143 0.0086 下 并 米 六 
草地 Grassland 31.56 0.2437 0.0264 0.0022 上 升 * 0.4703 0.0741 0.0067 在 外 水 水 
水 体 与 湿 } 
下 休 二 浊 地 3.77 0.2731 0.0281 0.0010 上 升 0.5315 0.0861 0.0096 此 升 玉米 
Water and wetland 
荒漠 Desert 13.5 0.0789 0.0043 -0.0002 下 降 六 冰 0.1543 0.0231 -0.0006 下 降 


* 为 趋势 通过 5% 显著 性 检验 , * * 为 趋势 通过 1% 显著 性 检验 


图 2 表明 ,农田 森林 和 草地 生态 系统 的 ANDVI 均 以 上 升 为 主 ,农田 生态 系统 上 升 像 元 比例 达 70.69% ， 
和 森林 和 草地 生态 系统 分 别 为 68.69% 和 63.49%。 人 然而, 上升 像 元 中 ,农田 、 和 森林 和 草地 生态 系统 分 别 仅 有 
39.86% .27.31% 和 27.19% 的 像 元 通过 了 5% 显著 性 检验 ,可 见 , 各 生态 系统 ANDVI 近 30 年 的 上 升 态势 并 非 稳 
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定 的 线性 趋势 ,而 是 在 上 升 过 程 中 伴随 着 较 大 的 波动 。 近 30 年 ,荒漠 生态 系统 有 73.61% 的 像 元 呈 下 降 趋 势 ， 
其 中 有 5$4.47% 的 像 元 通过 了 5% 显著 性 检验 ,荒漠 生态 系统 的 植被 存在 明显 的 退化 过 程 。 水 体 与 湿地 生态 系 
统 的 上 升 与 下 降 像 元 比例 几乎 相当 。 各 生态 系统 MNDVI 近 30 年 的 趋势 分 析 结 果 基 本 与 ANDVI 的 结果 一 
致 , 即 农 田 、 和 森林 和 草地 生态 系统 的 植被 呈 波 动 上 升 态 势 , 上 且 显 车 性 更 强 ,荒漠 生态 系统 呈现 较 明 显 的 退化 趋 
势 , 区 别 在 于 水 体 与 湿地 生态 系统 的 MNDVI 有 64.29% 的 像 元 表现 为 上 升 趋势 。 
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图 2 各 生态 系统 ANDVI 和 MNDVI 不 同 变化 趋势 的 面积 比例 
Fig.2 Area proportion of different trends on ANDVI and MNDVI in different terrestrial ecosystems 
FLE, 农田 生态 系统 farmland ecosystem; FOE, 森林 生态 系统 forest-ecosystem; GLE, 草地 生态 系统 grassland ecosystem; WWE, 水 体 与 湿地 生 
态 系统 water and wetland ecosystem; DE， 荒漠 生态 系统 desert ecosystem. | : 极 显 著 下 降 , ;显著 三 降 , 卫 ; 下 降 但 不 显著 , 玉 V: 上 升 但 不 显 
车,V : 显 车 上升 , YL: 极 显著 上 升 


2.3 不 同 陆 地 生态 系统 的 NDVI 波动 特征 及 残余 趋势 

对 近 30 年 中 国 各 陆地 生态 系统 的 NDVI 时 间 序 列 进行 集合 经 验 模 态 分 解 , 均 可 得 到 8 个 本 征 模 态 函数 
分 量 (IMF) 太 1 个 残 差 分 量 (RSE ) , 高 频 分 量 (IMEI 一 IMF3/4) 主要 为 模型 添加 的 日 电 声 和 植被 年 内 季 广 性 
生长 波动 信息 ,低频 分 量 (IMF4/5 一 IMF8 ) 为 年 际 \ 年 代 际 的 植被 波动 周期 ( 表 2)。 同 时 ,利用 日 噪声 的 集合 
扰动 ,检验 不 同 波动 周期 IMF 分 量 的 物理 意义 。 


表 2 1982 二 2012 年 不 同 陆地 生态 系统 NDVI 序列 各 分 量 的 平均 周期 
Table 2 Mean period of each intrinsic mode function of NDVI in different terrestrial ecosystems during 1982 一 2012 


生态 系统 类 型 Types of ecosystem 


本 征 模 态 函数 

Ta oc i 农田 条 林 草地 水 体 与 湿地 荒 沉 
Farmland Forest Grassland Water and wetland Desert 

4 2.0 1.9 / 7 % 

5 3 93 3.6 3.37 2.7 

6 ya 6.9** 6.8”* 6.8*” 5.3 
7 10.7* 11.7 10.3 10.9 9.6* 

8 35.3 36.7 22 .20 37.8 21.9 

# 表明 P<0.05, “表明 P<0.01;/ 表 明 该 IMF 不 具有 相应 尺度 的 周期 


由 于 大 气 与 海洋 环流 的 脉动 和 人 类 活动 间接 作用 于 陆地 生态 系统 ,导致 植被 活动 可 能 存在 周期 性 特征 。 
各 生态 系统 低频 分 量 的 平均 周期 ( 表 2) 和 显著 性 检验 (图 3) 表 明 , 农 田 生 态 系统 的 IMF 分 量 除 了 IMF8 外 , 均 
通过 了 5% 显 车 性 检验 ,可 见 农田 生态 系统 的 NDVI 具有 从 年 际 到 年 代 际 的 多 重 波动 周期 ,但 随 春 波动 周期 的 
增长 ,其 显著 性 逐渐 下 降 ,农业 生产 对 市 场 规律 的 适应 和 长 周期 的 气候 波动 都 可 能 是 影响 因素 。 和 森林 生态 系 
统 1.9a、5.3a 和 6.9a 的 波动 周期 通过 了 1% 显著 性 检验 ,但 长 周期 的 年 代 际 波动 显著 性 很 低 , 可 能 由 于 森林 生 
态 系统 的 群落 环境 相对 比较 优越 , 误 长 的 群落 演 蔡 过 程 受 气候 波动 影响 较 小 。 由 于 草地 生态 系统 包含 了 北方 
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温带 草地 .青藏 高 原 高 寒 草 地 及 南方 和 东部 次 生 草 地 ， 0 

区 域 气候 条 件 差异 大 ,植被 活动 的 气候 波动 痕迹 不 明 训 

显 ,该 生态 系统 仅 6.8a 和 22.2a 的 波动 周期 通过 了 5% i 

的 显著 性 检验 。 水 体 与 湿地 生态 系统 的 植被 活动 更 易 《| 。 

受到 气候 波动 的 影响 ,主要 表现 出 两 个 短 周期 的 波动 。 I 人 
荒漠 生态 系统 的 NDVI 可 能 包含 了 大 量 非 植被 的 背景 “| ---- 95% 显 著 性 
言 息 , 仅 9.6a 周期 通过 了 5% 显著 性 检验 ,可 以 认为 荒 

漠 生 态 系统 的 周期 特性 不 具有 明确 的 物理 意义 。 上 述 二 

分 析 表 明 ,以 上 5 类 生态 系统 主要 存在 1.9 一 7.6a 的 年 图 3 不 同 陆地 生态 系统 低频 分 量 显著 性 检验 


际 波 动 周期 特 征 和 9.6 一 11.7a 7 279 的 年 代 际 波 动 周 Fig.3  _ Significance test of low frequency component of NDVI in 
期 特征 ;其 中 年 际 波动 周期 特征 普遍 比较 显著 ,年 代 际 erent terrestrial ecosystems 
波动 周期 特征 显著 性 相对 较 差 。 由 此 推论 ,各 生态 系统 证 pyrion 
植被 活动 的 周期 性 变化 可 能 与 海洋 和 大 气 环流 的 脉动 有 一 定 联 系 。 

残 差 分 量 RSE 反映 NDVI 序列 的 内 在 发 展 趋 热 ,图 4 可 以 看 出 , 近 30 年 来 ,虽然 农田 、 森 林 .草地 和 水 体 
与 湿地 生态 系统 的 NDVI 呈 总 体 上 升 趋势 ,荒漠 生态 系统 的 NDVI 总 体 下降, 但 内 在 包含 了 “上 升 一 下 降 ” 两 
个 阶段 ,由 之 前 的 明显 上 升 趋 势 转变 为 不 同 程度 的 下 降 趋 势 , 这 与 传统 线性 趋势 显然 不 同 “” 。 暧 干 化 趋势 
下 ,尤其 是 生长 季 降 水 减少 引起 的 干旱 胁迫 是 植被 活动 发 生 转 折 的 重要 原因 ，””。 农田 森林 、 草 地、 水 体 与 
湿地 和 荡 漠 生态 系统 NDVI 趋势 转折 的 节点 有 所 差异 ,主要 集中 在 20 世纪 .90 年 代 和 本 世纪 初 ,植被 活动 转 
折 期 的 差异 在 相关 文献 “有 所 印证 ,这 可 能 是 气候 变化 响应 的 差异 和 生态 系统 分 类 精度 共同 作用 的 
结果 。 
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图 4 各 陆地 生态 系统 NDVI 残 差分 量 与 线性 趋势 


Fig.4 Linear trends and residual standard error on NDVI of every terrestrial ecosystem 
2.4 不 同 生态 系统 的 NDVI 空间 异 质 特征 分 析 
对 于 单一 生态 系统 的 NDVI 图 像 , 各 像 元 NDVI 之 间 的 标准 差 越 大 , 则 表明 该 生态 系统 植被 长 势 的 空间 异 
质 性 越 强 ,植被 活动 的 地 域 差 异 越 大 。 从 表 3 看 出 ,各 生态 系统 空间 异 质 性 也 存在 2.1 一 7.1a 的 年 际 波动 周期 
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特征 和 10.7a .30.2a 的 年 代 际 波动 周期 特征 ,但 各 生态 系统 空间 异 质 性 的 年 代 际 周期 波动 普遍 不 具有 物理 意 
义 ,主要 表现 为 2.1 一 7.1a 的 年 际 周期 节律 ,这 可 能 反映 了 短 周 期 气候 波动 对 生态 系统 植被 活动 的 影响 。 
表 3 1982 一 2012 年 中 国 陆地 生态 系统 NDVI 标 准 差 序列 各 分 量 的 平均 周期 


Table 3 Mean period of each intrinsic mode function on spatial heterogeneity of Chinese different ecosystems during 1982 一 2012 


生态 系统 类 型 Types of ecosystem 


本 征 模 态 函数 - 二 
Intrinsic mode function 农田 森林 单 地 水 体 与 湿地 抑 次 
Farmland Forest Grassland Water and wetland Deseff 

4 从 2 pl 2 A 
5 2.1 3.0 S90 四 人 3.6 
6 Sh Se ya a 3 

10.7 加 10.7 10.7 11.9 
8 16.7 29.8 32.3 21.9 302 一 


* 表明 P<0.05, ** 表明 P<0.01;/ 表 明 该 IMF 不 具有 相近 尺度 的 周期 


标准 差 序列 的 残 差 分 量 (RSE ) 表 明 ( 图 5) , 近 30 年 中 国 各 陆地 生态 系统 植被 长 势 的 空间 异 质 性 总 体 增 
强 , 植 被 活动 的 区 域 差 异 越 来 越 大 。 荒 漠 生 态 系统 NDVI 的 空间 异 质 性 持续 增强 , 直到 本 世纪 后 才 逐 渐 减 弱 ， 
这 与 亚洲 沙漠 边缘 植被 活动 变化 的 趋势 比较 一 致 ” ,沙漠 边缘 有 效 水 分 的 变化 直接 影响 到 植被 活动 ,进而 反 
映 到 空间 异 质 性 的 变化 趋势 。 和 森林 生态 系统 植被 长 势 的 空间 异 质 性 在 20 世纪 90 年 代 触 底 后 持续 增强 , 农 
田 .草地 和 水 体 与 湿地 生态 系统 的 空间 异 质 性 一 直 处 于 增强 趋势 区域 差 蜡 加 大 ,Peng 等 指出, 近 30 年 中 
国 南 北方 的 植被 长 势 存在 明显 差异 , 暖 干 化 趋势 对 中 国 北方 植被 长 势 的 季节 性 抑制 可 能 是 中 国 陆 地 生态 系 
统 植被 活动 空间 寞 奈 性 增强 的 重要 原因 。 


0219 -ee FLE <===“=FOE 二 二 人 [一 WWE 
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0.183 0.070 


残 差分 量 Residual standard error 


残 差分 量 Residual standard error 


0.165 0.065 
1982 1987 1992 1997 2002 2007 2012 
年 份 Year 
图 5 陆地 生态 系统 NDVI 标准 差 序列 的 残 差 分 量 


Fig.S Linear trends and residual standard error on spatial heterogeneity of every terrestrial ecosystem 
FLE, 农田 生态 系统 farmland ecosystem; FOE, 森林 生态 系统 forest ecosystem; GLE, 草地 生态 系统 grassland ecosystem; WWE, 水 体 与 湿地 生 
态 系 统 water and wetland ecosystem; DE, 蕊 漠 生 态 系统 desert ecosystem. 图 中 右 侧 纵 轴 仅 代表 总 漠 生 态 系 统 DE , 左 侧 纵 轴 代表 其 他 4 类 生态 
系统 


2.S、NDYVI 变化 周期 与 气候 事件 关系 分 析 

前 述 分 解 的 周期 与 ENSO 事件 2 一 7a 的 变化 周期 太阳 活动 有 关 的 黑子 周期 (11 一 12a) 海尔 周期 (22a) 、 
准 35a 的 布 虽 交纳 周期 在 时 间 上 存在 某 种 契合 ,可 能 表明 不 同 陆地 生态 系统 植被 活动 与 引起 海洋 和 大 气 环流 
脉动 的 ENSO 事件 ,太阳 活动 和 海 气相 互 作 用 等 存在 一 定 联系 。Propastin “和 Nagai “等 指出 陆 表 植 被 活动 
与 ENSO 事件 显著 相关 ,Li 等 利用 GIMMS 数据 研究 发 现 中 国 各 陆地 生态 系统 有 32.34% 的 区 域 的 植被 对 
ENSO 事件 响应 敏感 。 本 研究 将 MEI 指数 进行 EEMD 分 解 ,并 与 各 主要 生态 系统 NDVI 具有 相近 波动 周期 
的 IMF 分 量 进行 相关 分 析 ,发 现 ENSO 事件 存在 2.1a .5.3a 和 7.1a 的 波动 周期 ,生态 系统 NDVI 的 波动 周期 与 
ENSO 事件 存在 不 同 程度 的 相关 性 ( 表 4)。 农 事 生 产 通过 有 目的 的 水 分 调节 适应 农时 ,明显 掩盖 了 自然 气候 
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系统 对 农田 生态 系统 的 周期 性 影响 , 仅 IMF6 分 量 与 ENSO 事件 的 7.1a 周期 极 显 著 相关 ;森林 生态 系统 受 人 
类 活动 干扰 较 少 ,IMF4 JIMF5 和 IMF6 分 量 分 别 与 ENSO 事件 的 2.1a、5.3a 和 7.1a 周期 极 显著 相关 ,有 着 明显 
的 气候 短 周 期 脉动 痕迹 ;草地 生态 系统 的 周期 性 特征 不 明显 ,震荡 周期 与 ENSO 事件 无 关 ; 水 体 与 湿地 生态 系 
统 的 IMF6 分 量 与 ENSO 事件 7.1a 周期 极 显著 相关 ,湿地 植被 活动 存在 明显 振荡 周期 ;荒漠 生态 系统 植被 信 
息 不 足 ,波动 周期 与 ENSO 事件 无 显著 相关 性 。 

此 外 ,太阳 活动 与 全 球 气候 变化 关系 密切 “ ,尽管 太阳 活动 峰 、 谷 值 年 的 太阳 总 辐射 变化 极 小 ,但 施加 大 
气 系统 的 影响 可 以 通过 热带 地 区 海 - 气 耦合 作用 得 到 强化 ,Hadley/ Walker 环流 随 着 类 La Nifa 事件 的 出 现 而 
减弱 ” ,这 对 东亚 地 区 大 气 环流 异常 和 中 国 南北 地 区 反 向 降水 变化 影响 显著 “ ,而 降水 .气温 和 辐射 等 气候 
变量 则 通过 光合 作用 等 生理 过 程 影响 植被 的 生长 状况 和 群落 演 奉 。 太 阳 活 动 11a 周期 与 不 同 生 态 系 统 的 
NDVI 区 化 存在 不 同 程度 的 极 显著 负 相 关头 系 ( 表 4) ,这 可 能 与 太阳 活动 峰值 年 我 国 南方 和 东北 夏季 降水 偶 
少 有 关 ,而 生长 季 尤 其 是 夏季 降水 对 植被 生长 极为 重要 ,显然 植被 活动 11a 左右 的 震荡 周期 含有 明显 的 太阳 
活动 信息 。 上 述 分 析 表 明 ,气候 系统 的 脉动 会 引起 陆地 生态 系统 的 NDVI 波动 ,海洋 与 大 气 环 流 的 脉动 与 各 
生态 系统 植被 活动 的 周期 性 节律 有 着 一 定 关联 ,ENSO 事件 和 太阳 活动 是 推动 植被 活动 周期 性 振荡 的 因素 
< 


表 4 各 生态 系统 与 ENSO 太阳 黑子 相应 周期 分 量 相关 分 析 
Table 4 Correlations of the Intrinsic Mode Function between NDVI of each ecosystem and ENSO and sunspot 


气候 事件 周期 农田 森林 草地 水 体 与 湿地 荒漠 
Climate events Periods/a Farmland Forest Grassland Water and wetland Desert 
ENSO 2.1 -0.0888 -0:1634 全 于 = 一 
5.3 = -0.1306、* 三 = 0.0167 
7.1 -0.3147** -0.5455 ** -0.0701 -0.2529 
太阳 黑子 Sunspot 11 -0.2837 3 一 0.4439 ** 一 0.4038 ** —0.1761 ** —0.3180** 


ENSO ,厄尔尼诺 -南方 涛 动 El Niio-Southern Oscillation; * 显著 相关 (P<0.05), *##* 极 显著 相关 ( P<0.01) ;一 表明 该 生态 系统 不 具有 相近 尺度 
的 周期 


3 结论 


基于 1982 一 2012 年 的 GIMMS NDVl3g 和 中 国 陆 地 生态 系统 分 布 数据 ,利用 一 元 线性 回归 、 集 合 经 验 模 态 
分 解 和 相关 分 析 等 手段 ;研究 了 近 30 年 中 国 各 陆地 生态 系统 NDVI 的 时 空 变化 特征 ,初步 分 析 了 NDVI 变化 
与 气候 事件 的 关系 ,研究 取得 如 下 几 点 结论 : 

(1) 我 国 近 30 年 来 NBVT 主 要 以 上 升 为 主 ,ANDVI 的 上 升幅 度 为 0.0029/10a( P<0.05) ,MNDVI 的 上 升幅 
度 为 0.0076A10a( P<0.01)。 植 被 活动 显著 增强 的 区 域 主要 是 分 布 在 东部 季风 区 的 农田 和 和 森林 生态 系统 , 显 
著 下 降 的 区 域 主要 是 分 布 于 西北 的 荒漠 生态 系统 和 东北 森林 区 。 

(2) 从 线性 特征 来 看 ,尽管 各 ANDVI 和 MNDVI 变化 趋势 的 显著 性 有 所 差异 ,但 农田 、 和 森林 .草地 和 水 体 
与 湿地 生态 系统 的 NDVI 总 体 呈 现 非 稳定 的 上 升 趋势 ,上 升 过 程 中 伴随 着 较 大 波动 , 苇 漠 生态 系统 的 NDVI 呈 
下 降 趋势 ,植被 退化 显著 。 

(3)5 类 生态 系统 的 植被 活动 主要 存在 1.9 一 7.6a 的 年 际 周期 波动 特征 和 10.7a、22.2a 的 年 代 际 周期 波动 
特征 ,其 中 1.9 一 7.6a 的 年 际 周期 特征 普遍 比较 显著 ,年 代 际 周期 的 显著 性 相对 较 差 。 与 线性 趋势 不 同 ,各 生 
态 系 统 植 被 活动 的 残 差 趋势 包含 “< 上升 一 下 降 ” 两 个 阶段 ,并 相继 于 20 世纪 90 年 代 到 本 世纪 初 发 生 转 折 。 

(4) 近 30 年 来 ,我 国 各 生态 系统 植被 活动 的 空间 异 质 性 存在 2.1 一 7.1a 的 年 际 周期 市 律 ; 随 着 NDVI 的 上 
升 ,农田 和 森林、 草地 和 水 体 与 湿地 生态 系统 的 NDVI 空间 异 质 性 持续 增强 ,生态 系统 内 部 植被 活动 的 区 域 差 
异 逐 渐 加 大 , 苇 漠 生态 系统 植被 活动 的 空间 异 质 性 增强 的 趋势 于 本 世纪 初 后 逐渐 减弱 。 

(5) 气 候 系 统 的 脉动 会 引起 陆地 生态 系统 的 NDVI 波动 ,海洋 与 大 气 环流 的 短 周期 脉动 与 各 生态 系统 植 
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被 活动 的 周期 性 节律 有 一 定 关 联 ,ENSO 事件 和 太阳 活动 可 能 是 推动 植 彼 活动 周期 性 振荡 的 因 系 之 一 。 

太阳 活动 ,地球 轨 站 变化 ,造山 运动 和 人 类 活动 等 都 是 影响 地 球 气候 变化 的 重要 原因 ,而 太阳 活动 和 地 球 
轨 过 变化 表现 的 周期 性 节律 可 能 是 植被 活动 周期 性 振荡 的 重要 来 源 。 但 太阳 活动 和 地 球 轨 过 变化 与 气候 系 
统 的 联系 集中 在 相似 周期 的 比较 层面 ” ,其 影响 气候 系统 的 物理 机 制 并 不 十 分 明确 。 尽 管 本 研究 基于 集合 
经 验 模 态 分 解 初步 揭示 了 各 生态 系统 植被 活动 的 周期 性 振荡 特征 ,丰富 了 植被 活动 啊 应 气候 变化 的 认识 ,但 
该 方法 计算 的 平均 波动 周期 ”与 传 里 叶 谱 密度 的 结果 有 所 差异 , 且 生 态 系 统 内 部 植被 类 型 对 气候 变化 的 响 
应 差异 也 对 周期 的 显著 性 检验 存在 干扰 ,这 对 深入 解释 植被 活动 的 周期 特征 带 来 困难 。 近 30 年 来 ,中 国 大 范 
围 实施 了 退耕 还 林 .还 草 等 人 工 生态 治 理工 程 ,同时 还 经 历 了 快速 的 城市 化 ,人 为 因 系 引起 的 土地 利用 类 型 转 
变 也 对 各 生态 系统 的 植被 活动 有 着 深刻 的 影响 。 因 此 ,在 气候 变化 背景 下 和 剥离 生态 系统 转换 区 荆 植 被 活动 的 
变化 , 既 有 助 于 准确 量化 人 为 因素 对 植被 活动 的 影响 ,也 有 利于 解释 气候 变化 影响 植被 活动 的 机 制 。 
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